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We study level number variance in a two-dimensional random matrix model characterized by a 
power-law decay of the matrix elements. The amplitude of the decay is controlled by the parameter 
b. We find analytically that at small values of 6 the level number variance behaves linearly, with the 
compressibility \ between and 1, which is typical for critical systems. For large values of 6, we 
derive that \ — 0, as one would normally expect in the metallic phase. Using numerical simulations 
we determine the critical value of b at which the transition between these two phases occurs. 
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I. INTRODUCTION 

Random matrix models whose eigenstates exhibit a 
transition from extended to localized states provide an 
efficient tool for studying the Anderson metal-insulator 
transition [11-0] • Their main advantage, compared to the 
original Anderson model [8|, is that the transition oc- 
curs not just at a single point in the parameter space, 
but rather on a critical line described by the variation of 
an additional parameter, such as the band width of the 
power-law banded random matrix model Q. The mod- 
els are accessible to perturbative treatment when this 
parameter is either large or small [§-[l3j]. 

So far, most of the critical random matrix models 
which have been studied intensively are one-dimensional. 
Their Hamiltonians describe random hopping of a par- 
ticle on a one-dimensional lattice in a random on-site 
potential. The existence of the Anderson transition in 
such one-dimensional systems is related to the long-range 
nature of the hopping amplitudes. In our recent work 
[l4| , we studied the scaling of the moments of the eigen- 
states in a two-dimensional generalization of the power- 
law banded random matrix model [HI, [l|| . In this en- 
semble, the matrix elements of the Hamiltonian H ma 
are complex independent Gaussian random variables, 
whose mean values are equal to zero and whose vari- 
ances are determined by the distance between sites of a 
two-dimensional lattice: 



<|ff m n| 2 > = T- 
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where m = (m x ,m y ), n = (n x ,n y ), 1 < m a ,n a < L 
are two-dimensional vectors representing two sites on a 
two-dimensional square lattice of size L x L, and b is a 
parameter of the model. Thus the Hamiltonian is repre- 
sented by a random L 2 x L? matrix. 

One natural way to convince oneself that this system is 
critical is to study how the moments of its eigenfunctions 
V'n(r) scale with the system size L. Namely, we define 



where d is the dimensionality of the space and the av- 
eraging is performed over the ensemble as well as over a 
small energy window. Trivial exponent values d q = and 
d q = d signify localized and extended states, respectively. 
For critical states, < d q < d and what is more, d q gener- 
ally depends on q, indicating that the eigenfunctions are 
multifractal. The above scaling of the moments of the 
eigenfunctions was obtained analytically and confirmed 
by numerical simulations in various critical random ma- 
trix ensembles 0- In particular, it was found that in 
the power-law banded random matrix model d q -C 1 
when the bandwidth b <C 1 (strong multifractality) and 
d-(i 9 «l when b 3> 1 (weak multifractality) [!, [9| . 

One of the surprising features of the model (HJ revealed 
in [14j was the absence of a pure power-law scaling of 
the moments of the eigenfunctions at b ^> 1. Instead of 
Eq. @ , the scaling was surmised to be 
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where the exponents v q play the role of the anomalous 
fractal dimensions d — d q and can be calculated pertur- 
batively in a similar way In the same time, the 

standard power-law scaling ([2]) with d q -C 1 was found at 
b <C 1, in full analogy with the one-dimensional version 
of the model. 

The existence of two different scaling laws at large and 
small values of b suggests that there should be a critical 
value of b = 6 cr separating these two regimes. One of the 
aims of the present paper is to show that such a critical 
value of b does exist. While the model can be treated 
analytically at b <C 1 or b S> 1, the existence of the tran- 
sition between two regimes can be investigated only with 
the help of numerical simulations. These are more effi- 
cient for spectral properties rather than for the statistics 
of the eigenvectors. For this reason, in the present work 
we focus on studying the spectral compressibility. 

The spectral compressibility x is defined by the asymp- 
totic behavior of the level number variance: 

(5n 2 (E)) = (n 2 (E)) - (n(E)) 2 a X (n(E)), (n(E)) » 1, 

(4) 
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FIG. 1: Level number variance as a function of the average 
level number, for different values of 6 < I and L — 64. 



where n{E) is the number of eigenvalues in a spectral 
window of the width E. It is well known that x P iavs the 
role of a critical exponent: x = in the metallic phase, 
1 in the localized phase and < x < 1 at criticality 



Our main results concerning the behavior of the level 
number variance and the compressibility in the two- 
dimensional model (JlJ can be formulated as follows. The 
asymptotic behavior of (Sn 2 (E)) is described by Eq.Q 
for all values of b such that < b < b cr . The spectral 
compressibility is a monotonically decaying function 
of b with x(0) = 1 and x(^cr) = 0. The latter equation 
is used as the definition of 6 cr . At small values of b, 
can be calculated perturbatively: 



X(b) = 1 



V2 



o(b 4 ). 



(5) 



For b > 6 cr , spectral compressibility is equal to zero. 
The level number variance in this case has the same 
asymptotic behavior as in the Wigner- Dyson random ma- 
trix theory: 

(Sn 2 {E)) = i(ln(27r(n))+7 + l), (n) > 1, (6) 

where 7 is the Euler constant. Our analysis, based on 
numerical simulations, shows that the transition between 
these two regimes occurs at 6 cr = 5.2 ± 0.2. 

The remainder of this paper has the following struc- 
ture. In Section [IT] we present a derivation of Eq.© at 
b <C 1 and confirm the obtained result by numerical sim- 
ulations. The opposite case of b 3> 1 and the transition 
between the two phases are considered in Section IIIII 
We summarize our results in Section IIVI Finally, the 
Appendix contains an alternative derivation of Eq.©, 
easily to extended to the orthogonal symmetry class. 



FIG. 2: Spectral compressibility as a function of b, obtained 
from the slopes of the lines in Fig. [T] The data is presented 
for three different system sizes. The error bars are smaller 
than the symbol sizes. 



II. SPECTRAL COMPRESSIBILITY AT 6 < I 

A typical random matrix from the ensemble (|TJ) has 
the diagonal elements of order one and the off-diagonal 
elements of order b. Thus, the off-diagonal elements are 
parametrically smaller than the diagonal ones provided 
that b <C 1. This allows developing a perturbation expan- 
sion of the eigenfunction moments, spectral compressibil- 
ity or any other quantity of interest in a power series of 
b. In Ref.0, [l?} this method was applied to the calcu- 
lation of spectral compressibility in the one-dimensional 
power-law banded random matrix model. Below we use 
the same approach in the two-dimensional model ([I} . 

The zeroth order term in the expansion corresponds 
to a pure diagonal random matrix, which has completely 
localized eigenvectors and uncorrelated eigenvalues. As 
a result, x = 1 in this case. The first order correction to 
this trivial result can be found from the relation between 
X and the form factor K^(t,N)^: 



X = l + 



lim lim 

t— >oo N— >a 



(7) 



where N is the matrix size, which is equal to L 2 in our 
case, and the form factor is calculated in the lowest order 
of the perturbation theory. The general expression for 
the latter reads [TJ-E 7 ! 



2^ 
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s(|m-n|) = -(\H mn \ 2 )t 2 . 



(8) 



(9) 



Since we are interested in the limit N — > 00 , the sum over 
m in Eq.© can be replaced by a two-dimensional inte- 
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gral, which can be transformed to the polar coordinates 
upon substitution of (|-ff m n| 2 ) given in ([XJ) : 



lim K (1) (t,N) = -2tt 3/2 

N— ¥00 



/; d w (10) 



This integral can be evaluated in the limit t 
changing the variable from r to s — bH 2 /2r 4 : 



00 by 
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V~s~ s/2' 
(11) 

leading to the formula announced in Eq.([5]). The same 
result can be obtained directly by calculating \ for a ran- 
dom 2x2 matrix, as shown in the Appendix. The cor- 
responding expression for \ in the orthogonal symmetry 
class (real symmetric matrices) is also given there. 

In order to test this prediction we performed numerical 
simulations with random matrices generated according to 
Eq. (JTJ) . The spectra of these matrices were obtained by 
standard diagonalization subroutines [20], and the num- 
ber of realizations for each b was 20000, 2000 and 500 for 
L = 64, 128 and 180, respectively. 

The numerical findings for b < 1 are summarized in 
Fig. [TJ and Fig. [2] The first one shows the level number 
variance (Sn 2 (E)) as a function of (n(E)) for different 
values of b. The linear behavior suggested by Eq.Q is 
evident for b < 1. Extracting the slopes of the straight 
lines we obtain the numerical values for x(b), which are 
presented in Fig. [5] They are in a good agreement with 
our analytical prediction ((SJ, which is valid only when 
the absolute value of the perturbative correction ([TT]) is 
much smaller than one, i.e. b 2 <C V2/n 2 ps 0.14. 



III. LEVEL NUMBER VARIANCE AT b > 1 
AND THE TRANSITION POINT 

To study the limit b 3> 1 we exploit the fact that in 
this case the model can be mapped onto a non-linear a- 
model flij . In the a- model description the key ingredient 
characterizing a particular model is the propagator. In 
Ref.[3] it was found that the propagator of the two- 
dimensional model (TTJ) is given by 



n(fc) 



1 



2 k 2 \nbk 



(12) 



in the momentum space. Its eigenvalues determine vari- 
ous spectral properties of the model [2l], [22[ • In partic- 
ular, they allow us to calculate the two-level correlation 
function 



R{u ) = ^Re 

k 



1 



^-k 2 \xxbk-iu 



(13) 



where A is the mean level spacing, v = 1/AL 2 is the 
density of states and the sum runs over the discrete mo- 
menta k = (2irm x /L 1 2irm v /L), m a <= Z. The two-level 



correlation function allows us to find the derivative of 
(5n 2 (E)) [HI 



d(6r 



d(n) 



r{n) 

= lim / ds R(s), s = u>/A. 

L-yooJ_, n) 



(14) 



Substituting the expression (|13[) into this formula and 
integrating over s we obtain 



d(Sn 2 ) _ (to) 
d(n) ~ (8tt6 2 ) 2 



lim > 

r v^, Z — / 



1 



(»)' 



(15) 

Since the sum over m = (m Xy m y ), m a € Z converges 
and each term with m ^ (0, 0) tends to zero in the limit 
L — > 00, we conclude that only the zero mode k = has a 
non-vanishing contribution in the thermodynamic limit. 
It is well known that the contribution of the zero mode of 
the er-model reproduces the results of the Wigner-Dyson 
random matrix theory [23]. Thus we must expect that 
in this regime (Sn 2 (E)) is given by the standard random 
matrix theory result |24| : 

(Sn 2 {E)) = -l(ln(27r(ra))+7 + l), (n) > 1. (16) 

The absence of the linear in (to) term in this asymptotic 
law implies that \ = 0, as one usually finds in the metallic 
phase. We would like to point out that similar calcula- 
tions for the one-dimensional version of the model give 
a non-zero result for x S Si] ■ It is the logarithmic term 
in the propagator (|12j) . that leads to the vanishing of \- 
The same logarithmic term is responsible for the unusual 
scaling behavior of the moments of the eigenfunctions 

Fig. [3] shows the results of numerical simulations for 
b > 1. We observe a very good agreement between Eq. 
([T6"f and the numeric data at b = 5. One can also see how 
the behavior of (8n 2 ) changes gradually from linear for 
small b to logarithmic for large b. To quantify this change, 
as well as to determine the transition point between the 
two regimes, we assume that the most general functional 
form of (Sn 2 ) is given by 



(5n 2 ) = x(n) 



Oi 



ln(27r(n)) +a 2 , (to) > 1. (17) 



This formula naturally interpolates between Eq. Q and 
Eq. (|16p . Moreover we note that Eq.Q is never realized 
in its pure form: there is always a subleading logarithmic 
term. If the transition between the two regimes occurs 
at some finite value of b — 6 cr , one should expect that \ 
decays as a function of b and goes to zero at 6 cr . 

Using x, a± and 02 as fitting parameters, we were able 
to reproduce the behavior of (Sn 2 ) at intermediate values 
of b as shown in Fig- [3] Moreover, the best- fit parameters 
for large values of b are in agreement with Eq.(fT6"|). At 
6 = 5, for example, we find X = 0.00012 ± 0.00016, 01 = 
1.007±0.007 and a 2 = 0.1573±0.003, while the standard 
random matrix theory prediction (|16p corresponds to % = 
0, ax = 1 and a 2 = 0.1598. 
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FIG. 3: Level number variance as a function of the average 
level number, for different values of b > 1 and L — 64. Dots 
are numerical data. The solid line at b — 5 is Eq. (|16[) . The 
other solid lines are the results of fitting the data with Eq. 

G3- 



The values of Xi obtained in this way for different val- 
ues of b, are presented in Fig.|U One can see that, indeed, 
X monotonically decreases as a function of b and becomes 
zero within numerical accuracy for b > 6 cr . A close study 
of x in the interval 4.5 < b < 5.5 allowed us to determine 
the critical value 6 cr = 5.2 ± 0.2. 



IV. CONCLUSIONS 

The random matrix model described by Eq.([T]) was 
studied numerically in Ref. [HI, [l6[ . The results of those 
works suggest that the model is critical at all values of b. 
All analytical and numerical findings known for the one- 
dimensional counterpart of the model support the same 
expectation @, 0|- 

In this paper we show that the situation is actually 
more subtle. For 6< 1 we found the expected critical be- 
havior, characterized by a non-zero value of the spectral 
compressibility (JSJ) - Also, the scaling of the eigenstates 
is the standard power-law (0]) with non-trivial multifrac- 
tal dimensions d q . For b ^> 1, on the other hand, we 
obtained x — 0. This is normally a signature of the 
standard metallic phase corresponding to completely ex- 
tended states with d q = d. However, this is not the case 
in our model, where the moments of the eigenstates con- 
tain an additional anomalous part which scales as a power 
of the logarithm of system size ©. Therefore, the phase 
at 6 ^ 1 is not entirely metallic. We may say that it 
has some traces of the critical behavior. We found that 
the transition between the two phases - "critical" and 
"metallic-critical" - occurs at 6 cr = 5.2 ± 0.2. 




b 



FIG. 4: Spectral compressibility as a function of b, extracted 
from the data in Fig. [3] using Eq. (|17|l . The error bars are 
smaller than the symbol size. 

We believe that two important features of the model 
are responsible for the emergence of the two phases: the 
dimensionality d = 2 and the long-range nature of the 
Hamiltonian ([IJ. Recently, a similar behavior was pre- 
dicted for another two-dimensional long-range random 
Hamiltonian (2o| . 

Finally, we would like to point out that y given in 
Eq.© and di = \f2-K 2 b 2 , as calculated in Ref. |14j . satisfy 
the relation 

X+^ = l, (18) 

which was suggested recently in Ref. [26] and verified for 
various one-dimensional random matrix models. 

V. ACKNOWLEDGMENTS 

We are grateful to Yan Fyodorov, Vladimir Kravtsov 
and Alexander Mirlin for useful discussions. IR 
and AO acknowledge support from the Engineering 
and Physical Sciences Research Council [grant num- 
ber EP/G055769/1]. E.C. thanks the FEDER and the 
Spanish DGI for financial support through Project No. 
FIS2010-16430. 



VI. APPENDIX. CALCULATION OF 
SPECTRAL COMPRESSIBILITY AT b < 1. 

In the limit of an almost diagonal Hamiltonian the 
problem can be reduced to a 2 x 2 Hamiltonian Q: 
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whose eigenvalues A± are the solutions of the character- 
istic equation 

(e 1 -X)(e 2 -X)-b 2 \h\ 2 = 0. (20) 

Both €1,2 are Gaussian with variance a 1 . The density- 
density correlator of this 2x2 matrix correlator near 
E = is R(w) = \ {5(u - X + )6(X-)) . We are inter- 
ested in the true density correlator of the original large 
random matrix. Although the correlators in these two 
cases are related, they are not the same: it turns out 
that we need to double the 2 x 2 R(ui). Indeed, in 
a very large matrix diagonal matrix (the limit 6 = 0) 
(p(oj)p(O)) = (p(w)) (p(0)) , whereas in the 2x2 matrix 
the r.h.s. contains a coefficient 1/2. Doubling the corre- 
lator, we have 

R ( u ) = ^2 (/ e-^^w - A+)<J(A_)dei<fe 2 ^ , 

where Qh is the average over the off-diagonal elements. 
Using the Eq.. (|2"U|) we can change the integration vari- 
able ei —> A_. One integral is then removed by <5(A_). 

b 2 \h\ 2 

The other S- function ensures H — — = to, which can 

be used in the exponential. Afterwards it becomes clear 
that the exponential can be neglected, as long as we are 
working in an energy window u <C <x. In this limit, 

m = (-^ww lH>2Wl ) h (22) 

The average density of states at E — is given by a sim- 
pler and similar calculation: (p(0)) = , . Knowing 
the density correlator in the 2x2 case, we can find the 
leading order of the density correlator in the original large 
system by replacing bh with the off-diagonal elements Hij 
and introducing the summation over Then we can 
substitute the connected correlator into Eq. (fT4| . As a 
result, 

X = f — Urn I 2 (n) — 

(n) — voo \ 



2(n)^/ / 27V 2 I , ~ \\ . . 

^E<^-^^I^I 2 1,^ |2 <^^>)- (23) 

ijtj v 

This result is general - valid for any dimensionality and 
for any symmetry class. Note that each term of the sum is 
an independent random variable and so can be averaged 
separately. Performing the averaging in the unitary case, 
we obtain x = 1 — hm^^^ / where 

(24) 

Here we used the standard notation for the variance 
(|-ffr,r'| 2 ) = a rr'- ^ n our system it is sufficient to 
set a r , r ' = | r _^ r /|2 ■ To simplify the calculation, we can 
imagine that our system occupies a very large disk of 
radius R (so that N = ttR 2 ), replace summation with 
integration over r and r' and subsequently change the 
integration variables to r — r' and r + r'. The integrals 
can be computed in polar coordinates. The first inte- 
gration depends weakly on the large upper limit and so 
we can set it to R. Afterwards the integrand does not 
depend on r + r', making the second integration trivial. 
The result is 

Taking the large (n) limit, 

X (U) = l-^ + 0{b% (26) 

The case of the orthogonal symmetry is treated in pre- 
cisely the same way, starting from Eq. (1231) . There we 
obtain 

X {0) = l-2irb 2 + 0{b 4 ). (27) 
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